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ABSTRACT 

Context. The MHD instabilities can generate complex field topologies even if the initial field configuration is a very simple one. 
Aims. We consider the stability properties of magnetic configurations containing a toroidal and an axial field. In this paper, we 
concentrate mainly on the behavior of non-axisymmetric perturbations in axisymmetric magnetic configurations. 
Methods. The stability is treated by a linear analysis of ideal MHD equations. 

Results. In the presence of an axial field, it is shown that the instability can occur for a wide range of the azimuthal wavenumber 
m, and its growth rate increases with increasing m. At given m, the growth rate is at its maximum for perturbations with the axial 
wave-vector that makes the Alfven frequency approximately vanishing. We argue that the instability of magnetic configurations in the 
ideal MHD can typically be dominated by perturbations with very short azimuthal and axial wavelengths. 
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1. Introduction 

A wide variety of MHD instabilities can occur in magnetized 
astrophysical bodies where they play an important role in the 
evolution and formation of various structures and in enhancing 
transport processes, among others. The onset of instabilities can 
be caused both by hydrodynamic motions (for instance, differ- 
ential rotation) or properties of the magnetic configuration. Even 
magnetic fields with a relatively simple topology (for example, 
a purely toroidal field) can be subject to instability. Magnetic 
fields generated by the dynamo action or stretched by hydrody- 
namic motions are topologically more complex and can cause 
this sort of instability as well. Which field strength and topology 
can sustain a stable magnetic configuration is still rather uncer- 
tain despite all the extensive work already done (see Borra et al. 
1982; Mestel 1999 for review). 

The simplest and best-studied magnetic configuration is 
most likely a purely toroidal one. This has been known since the 
paper by Tayler (1957), where stability properties of the toroidal 
field B v are determined by the parameter a = dlnB^/dln s 
where s is the cylindrical radius. The field is unstable to axisym- 
metric perturbations if a > 1 and to non-axisymmetric pertur- 
bations if a > -1/2. The growth time of instability is close to 
the time taken for an Alfven wave to travel around the star on 
a toroidal field line. Numerical modeling by Braithwaite (2006) 
confirms that the toroidal field with B v oc s or oc s 2 is unstable to 
the m - 1 mode (m is the azimuthal wave number) as predicted 
by Tayler (1957, 1973). However, even a purely toroidal field can 
be stable in the region where it decreases rapidly with s. A purely 
toroidal field cannot be stable through the whole star because the 
stability condition for axisymmetric modes (a < 1) is incompat- 
ible with the condition that the electric current in the z direction 
has no singularity at s — » 0, which implies a > 1. The stabil- 
ity of the toroidal field in rotating stars has been considered by 
Kitchatinov & Rudiger (2007), who argue that the magnetic in- 



stability is essentially three-dimensional and that the finite ther- 
mal conductivity creates a strong destabilizing effect. Terquem 
& Papaloizou (1996) and Papaloizou & Terquem (1997) consid- 
ered the stability of an accretion disk with the toroidal magnetic 
field and found that the disks containing a purely toroidal field 
are always unstable and calculated the spectra of unstable modes 
in the local approximation. 

The stability properties of purely poloidal magnetic fields 
are also well-studied. It has been understood since the papers 
by Wright (1973) and Markey & Tayler (1973, 1974) that the 
poloidal field is subject to dynamical instabilities in the neigh- 
borhood of points (or lines) where the poloidal field is vanish- 
ing (neutral points/lines). These authors recognized first that the 
magnetic field in the neighborhood of a neutral line resembles 
that of a toroidal, pinched discharge, which is known to be un- 
stable. The instability of a poloidal field is also rather fast: its 
growth time can reach a crossing time of few Alfven times (Van 
Assche et al. 1982; Braithwaite & Spruit 2006) that is very short, 
for example, compared to the time-scales of stellar evolution. 
However, the instability of a poloidal field can be suppressed by 
the addition of a toroidal field in the neighborhood of neutral 
points (Markey & Tayler 1973; Wright 1973). 

Conversely, the addition of even a relatively weak poloidal 
field alters the stability properties of the toroidal field substan- 
tially. For example, if the poloidal field is uniform and rel- 
atively weak, the instability condition of axisymmetric modes 
reads a > -1, at variance with the condition of instability for 
a purely toroidal field (see, e.g., Knobloch 1992; Dubrulle & 
Knobloch 1993), which predicts that an unstable toroidal field 
configuration has a > -1/2. Therefore, a weak poloidal field 
has a destabilizing effect. However, a strong enough poloidal 
field can suppress the instability of the toroidal field. It turns out 
that configurations containing comparable toroidal and poloidal 
fields are more stable than purely toroidal or purely poloidal 
ones (Prendergast 1956; Tayler 1980) and, generally, the possi- 



bility exists that there are configurations containing mixed fields, 
which have no instabilities arising on a dynamical time-scale. 
In his study of unstable magnetic configurations Tayler (1980) 
has not found any instability if the axial field B z > 0.3Z^ for 
instance, even though such configurations can be unstable for a 
wide range of the azimuthal wavenumber m if B z is weaker. With 
numerical simulations Braithwaite & Nordlund (2006) studied 
the stability of a random initial field in the stellar radiative zone 
and found that the stable magnetic configurations generally have 
the form of tori with comparable poloidal and toroidal field 
strengths. 

In this paper, we consider in detail the stability properties of 
magnetic configurations containing the toroidal and axial mag- 
netic fields with respect to non-axisymmetric perturbations. We 
show that the instability may occur for a wide range of the az- 
imuthal wavenumber m, and the growth rate is typically higher 
for higher m. Unstable modes with large m have a very short ver- 
tical lengthscale, so it can be hard to resolve them in numerical 
calculations. Depending on the profile B<p(s) and the ratio B z /B <p , 
the instability can occur in two regimes that have substantially 
different growth rates. 

The remainder of this paper is arranged as follows. In Sect. 
2, we derive the equation that governs the eigenfunctions and 
eigenvalues of the magnetic field. We describe the numerical 
procedure and present the results of calculations in Sect. 3. A 
brief discussion of the results is given in Sect. 4. 

2. Basic equations 

Let us consider the stability of an axisymmetric cylindrical mag- 
netic configuration in a high conductivity limit. We work in 
cylindrical coordinates (s, if, z) with the unit vectors (e s , e v , e z ). 
We assume that the azimuthal field depends on the cylindrical 
radius alone, B v = B^s), but the axial magnetic field B z is con- 
stant. 

In the incompressible limit, the MHD equations read 



dv Vp 1 

— + (v • V)v = — - + (V xB)xB, 

at p 4np 

V • v = 0, 

^-Vx(vxB) = 0, 
at 

V • B = 0. 



(1) 
(2) 
(3) 
(4) 



In the basic state, the gas is assumed to be in hydrostatic equi- 
librium, then 



Vp 1 

— = (VxB)xB. 

p 4np 



(5) 



In this paper, we consider the stability of non-axisymmetric per- 
turbations. Since the basic state is stationary and axisymmetric, 
the dependence of perturbations on t, <p, and z can be taken in 
the form exp (erf - ik z z - imip) where k z is the wave-vector in 
the axial direction and m the azimuthal wavenumber. In stellar 
conditions, such a local analysis in the z-direction applies if k z 
satisfies the inequality k z s > 1 (we assume that the lengthscale 
in the axial direction is ~ s in stars). Small perturbations will be 
indicated by subscript 1, while unperturbed quantities will have 
no subscript. Then, the linearized Eqs. (l)-(4) are 



Vpi 1 

on = — — + [(V x Bi) X B + (V x B) x Bi], 

p 4np 



(6) 



Vvi =0, 

crBi -(B V)vi +(n • V)B = 0, 



VB, =0. 



(7) 
(8) 
(9) 



Eliminating all variables in favor of the radial velocity perturba- 
tion vis, we obtain 



d_ 

ds 



idvi s vu 

- r (cr z + ^)|- 7 - + — 



ds 



- k 2 z {(T 2 + w\)V U 



where 



2 m(l+A) 
lqco B (l - a) ■ 



s 2 A 2 
ma>Az 



(l - -^jj(co Az + 2mco B ) 



s 2 A 2 



Ak 2 (jj 2 A u) 2 B 

M^+CO 2 ) 



Vis = 0, 



(10) 



1 / m \ 

coa= , [k z B z + —B lfi \, a> Az = 
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4np 
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4np 



co B -- 



s -\J4ttp 



_ dlnB v ^ _ 1 + m 2 
d\ns ' s 2 k 2 ' 



For axisymmetric perturbations (m = 0), Eq. (10) recovers 
Eq. (11) of the paper by Bonanno & Urpin (2007), if one as- 
sumes B z =const. The equations derived by Acheson (1973) and 
Knobloch (1992) for axisymmetric perturbations can also be re- 
covered from Eq. (10). With appropriate boundary conditions, 
Eq. (10) allows the eigenvalues cr to be determined. 

Our analysis does not include gravity, which is important in 
stars, so the stabilizing effect of stratification is neglected. In the 
case of magnetic instabilities, this is justified if the work done 
by a perturbation against gravity is less than the energy released 
from the magnetic field. The corresponding condition for stel- 
lar radiative zones has been considered by Spruit (1999) (see 
Eq. (44)) and reads in our notations as k z s > N/co B , where N is 
the buoyancy frequency. Assuming that the characteristic value 
of ./V in stars is ~ 10~ 3 s _1 and introducing the axial wavelength 
A z = 2n/k z , we obtain that the effect of stratification can be ne- 
glected for perturbations with 



A z < A c = 2 X 10 s B^p JJ 2 cm, 



(11) 



where B v3 = B^/10 3 G and p_ 4 = p/10~ 4 g/cm 3 . 

Magnetic perturbations can also be affected by the presence 
of dissipation, which is not considered in our analysis. Usually, 
magnetic diffusion is weak in stellar conditions, and its influence 
is unimportant if the dissipation rate that is ~ r/k 2 is small com- 



pared tow s when r\ 
condition yields 



10 3 cm 2 /s the magnetic diffusivity. This 



A z > A d = 3 X 10 5 B-^ 2 s\[ 2 p-_\ /4 cm, 



(12) 



where iii=5/10 u cm. Inequalities (11) and ( 1 2) set a range of 
axial wavelengths where our consideration is valid. 

3. Numerical results 

We assume that the dependence of the azimuthal magnetic field 
on s is given by 



Bui - Bm — 



(13) 



where B^o is the field strength at s = si. To calculate the growth 
rate of the instability, it is convenient to introduce dimensionless 
quantities 



B, 



x=—, q=k z si, r= , w B0 = - 

*1 ^bo si ^4np 



Then, Eq. (10) transforms into 

dv\ s vi s \ din A 



s=-^. (14) 



d I dv\ s vis 
dx\ dx x 



dx 



dx 



<t 1 + 



2q 2 x a - 1 
F 2 +f 2 



q 2 x 2 



qx 2 



(I -a)- 



q 2 x 2 
2mf 



m 2 + q 2 x 2 



Vis ~ 



4q 2 f 2 x 2(a-l) 
+ (r 2 +/ 2)2 ^ = °< 



where 



f = qs + mx a 1 , A = 



q 2 x 2 (T 2 +f 2 ) 
m 2 + q 2 x 2 



(15) 



(16) 



Solving the eigenvalue problem numerically, we assume that the 
inner and outer boundaries correspond to xi = 1 and X2 = 2, 
respectively, where we suppose Vi.,(xi) = v\ s (x2) — 0. Note, 
however, that the results are not qualitatively sensitive to the 
choice of boundaries and boundary conditions. We checked that 
the qualitative behavior of the growth rate does not change if 
the outer boundary X2 is set to higher values. As far as the inner 
boundary is concerned, our approach does not allow the eigen- 
values to be calculated in the case si = 0, but we verified that 
low but finite values of s\ did not lead to qualitative difference 
in the behavior of the growth rate, and therefore we believe that 
our conclusions can be extended to the case si = 0. 

Equation (15), together with the given boundary conditions 
at the extrema, is a two-point boundary value problem that can 
be solved by using the "shooting" method (Press et al., 1992). 
To solve Eq. (15), we used a fifth-order Runge-Kutta integra- 
tor embedded in a globally convergent Newton-Rawson iterator. 
We have checked that the eigenvalue was always the fundamen- 
tal one, as the corresponding eigenfunction had no zero except 
that at the boundaries. To test the numerical procedure, we con- 
sider first the case of a purely toroidal field, e — 0. Stability 
properties of such field have been well studied since the paper 
by Tayler (1973), who argued that the dominant unstable mode 
is the one with m = 1 . The instability of the toroidal field oc- 
curs if d(sB 2 )/ds > or, in other words, if B v increases with the 
cylindrical radius or decreases not faster than s^ 1 ^ 2 . In Fig. 1, we 
plot the growth rate as a function of the vertical wave-vector for 
B, = and different values of a. In agreement with the conclu- 
sion by Tayler (1973), a = -1/2 is indeed the critical value dis- 
tinguishing between stable (a < -1/2) and unstable (a > -1/2) 
toroidal magnetic fields. The growth rate of instability can be 
rather high, cr ~ a> B o- Note that the growth rate increases with 
increasing q, and perturbations with short axial wavelengths turn 
out to grow the most rapidly. 

In Fig. 2, we plot the dependence of T 2 on q for the same 
m = 1 mode but in the presence of a relatively weak axial field, 
B z = 0.1 B^o- The addition of even very weak B z changes the sta- 
bility properties qualitatively even though the energy contained 
in the axial field is very low compared to that of the toroidal 
field (~ 1%). We stress that the presence of an axial field breaks 
symmetry q — > -q. The physical reason for this is fairly sim- 
ple. The Lorentz force plays a crucial role in the behavior of 




Fig. 1. The dependence of F 2 on q for B z = and m = 1. 
Numbers near the curves indicate the values of parameter a. 
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Fig. 2. The dependence of T 2 on q for s = B z JB ip o = 0.1 and 
m-\. Numbers near the curves correspond to values of a. 



perturbations, and this force contains a component that is pro- 
portional to the cross production of a current flowing in the ba- 
sic state j oc V x B and the magnetic field of perturbations (see 
the last term on the left hand side of Eq. (6)). If B z + 0, then 
magnetic perturbations are determined partly by the axial gradi- 
ent of velocity perturbations since dB\/dt = (BV)vi, and this 
contribution has a different sign for positive and negative k z . 
Therefore, the stability properties turn out to be dependent on the 




Fig. 3. The dependence of F 2 on q for s = BJB^ = 0.1, 
a = 0, and different values of m. Solid, long-dashed, dotted, 
short-dashed, and dot-dashed lines are for m = 0, 1,2, 3,4, re- 
spectively 

direction of an axial wavevector. However, Eq. (15) still contains 
some degeneracy because the replacements (m, q) — > (— m, -q) 
or (m, e) — » (— m, — e) do not change its shape. The instability oc- 
curs only for a restricted range of negative q, -(20-30) < q < -2 
depending on the value of a, and do not appear for positive q. 
The growth rate has two clear maxima with the higher maxi- 
mum corresponding to q 1 /s. By the order of magnitude, the 

axial wave-vector of the most rapidly growing perturbation can 
be estimated from the condition 

1 / m \ 

<o A = -—(k z B z +-BJ*0. (17) 

Indeed, this equation in dimensionless units reads a> A oc qs + 
mx ff ~' « 0. Since x ~ 1 in our calculations, the condition aj A = 
corresponds to 

q ~ -m/s. (18) 

Therefore, we have q ~ -1/e for the m — 1 mode. The most 
rapidly growing modes turns out to be highly anisotropic if the 
axial field is weak compared to the toroidal one. Their axial 
wavelength A, = 2n/k z ~ 2nss is much shorter than the radial 
and azimuthal lengthscale. The presence of an axial field shifts 
the threshold of unstable values of a: a purely toroidal field is 
unstable only if a > -1/2, but our calculations show that the 
profile of B<p with a = -1/2 is unstable. However, the growth 
rate in the case a < -1/2 is rather low and is close to to Az rather 
than <jj bo . The growth rate grows rapidly with increasing a: the 
maximum value of T 2 is of the order of unity for a = 1 and ~ 8 
for a = 2 (is not shown in Fig. 2). 

The dependence of F 2 on q in the presence of a weak axial 
field (e = 0.1) and for different values of m is shown in Fig. 3. 
In this figure, T 2 is calculated for a = 0, such that the magnetic 



Fig. 4. The same as in Fig. 3 but for B z /B v q = 0.5. 

configuration is unstable even in the absence of the axial field. It 
turns out that this profile of B^ is unstable if B z + 0, as well, but 
the properties of instability are substantially different. In contrast 
to the case of a purely toroidal field with a = where only the 
m — 1 mode can arise (Tayler 1973), the instability occurs for 
modes with a wide range of m including m — 0. For m + 0, the 
growth rate reaches its maximum approximately at q as -m/e, 
which is equivalent to the resonant condition oj a « 0. The max- 
imum growth rate increases with an increase in m. Therefore, 
the instability of such magnetic configurations is probably domi- 
nated by the modes with large m and extremely short axial wave- 
lengths A, as -2nes/m. This can cause difficulties in numerical 
simulations, because very high resolution in the axial and az- 
imuthal direction would be needed to resolve the most unstable 
modes. Moreover a nonlinear interaction of these modes can pro- 
duce axial lengthscales that are even shorter than those predicted 
by the linear theory. 

An increase in the axial field makes the magnetic configura- 
tion more stable. In Fig. 4, we plot the growth rate for s = 0.5 
and the same values of other parameter as in Fig. 3. The insta- 
bility grows more slowly in higher B z : an increase in B z by a 
factor 5 leads to the decrease in F 2 approximately by a factor 2. 
Nevertheless, the instability can still occur for this e, and it is 
still efficient because its growth rate is ~ 0.3ol>bq- Like the pre- 
vious case, the modes with a wide range of m can be unstable, 
and the growth rate increases with m. The axisymmetric mode 
(m = 0) is stable in this case (see Bonanno & Urpin 2007). The 
critical value s that suppresses the instability is clearly depen- 
dent on the geometry of the basic state. For example, by using 
energy considerations, Tayler (1980) has not found any instabil- 
ity for s > 0.3 in the configuration where the magnetic surfaces 
of the poloidal field are coaxial tori. On the other hand, for val- 
ues of s significantly less than 0.3, the author found instability 
for a wide range of values of m. 

To illustrate the dependence of F on the profile of the toroidal 
field, we show the growth rate as a function of q for a = 1 in 



9.0 [ 



1.2 



1.0 



\ ' \ 



7 \ 



0.6 



0.4 



0.2 



0.0 h^A 



i i 



A 



1 / 



'i ' ■ 

:i / i 



-60 -50 -40 -30 

q 



-20 



10 



8.5 



8.0 



7.5 



7.0 



6.5 



6.0 I 



/ \ 

i \ 



i ' i/ !'• 

i! ' :i ': 

II I j I/. V 

III , II 

I ■" . : l i 

I :, * 



I i i i i I ii i i III i ■ 



Fig. 5. The same as in Fig. 3 but for a = 1 . 
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Fig. 6. The same as in Fig. 3 but for a = 2. The dash-dot-dotted 
line represents m = 10. 



Fig. 5. This profile seems to be particularly interesting for astro- 
physical applications since the toroidal field is oc s in a neigh- 
borhood of the axis of symmetry. In the absence of the axial 
field, the profile with a = 1 corresponds to the toroidal field 
that is marginally stable to axisymmetric perturbations (m = 0). 
If B z + 0, the threshold of instability can change substantially 
as argued by Bonanno & Urpin (2007) and, indeed, perturba- 
tions with m = turn out to be unstable in this case. However, 
the non-axisymmetric perturbations grow faster and are likely to 
dominate the instability. The maximum growth rate is very high 
for them and comparable to the Alfven frequency for the toroidal 
field, (x>bq- Remarkably, the maximum growth rate changes rather 
slowly with m for all modes with m > 2. As in the previous cases, 
the maximum of F is located at q ~ -m/s. 

In Fig. 6, we show T 2 for the case of a rapid increase in B v 
with s, By oc s 2 or a = 2. The growth rate of instability is sub- 
stantially higher for such a and can reach the value ~ 2 - 3 
toroidal Alfven frequencies. A qualitative behavior of F remains 
same: the maximum growth rate is higher for modes with higher 
m, and these maxima correspond to very high values of the wave- 
vector q m/s. The maximum growth rate of the axisymmet- 
ric mode (m = 0) is comparable to that of non-axisymmetric 
ones (see Bonanno & Urpin 2007). Note that our calculations 
show some trend in the growth rate to reach saturation for large 
m. 

Figure 7 plots the growth rate for a = -0.5. In accordance 
with Tayler (1973), such a profile of B^ should be marginally 
stable if the axial field is vanishing. However, the presence of B z 
makes the toroidal field unstable with such a even if the axial 
field is relatively weak. The growth rate is not high in this case 
and is determined by the Alfven frequency for the axial field. 
As usual, the maximum growth rate is higher for higher m, and 
these maxima are reached for very large q, which corresponds to 
a short axial wavelength. A comparison between Figs. 6 and 7 
illustrates the difference well between two regimes of the insta- 
bility first noted by Bonanno & Urpin (2007). If a < a c where 
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Fig. 7. The same as in Fig. 3 but for a = -0.5. 



a c ~ 1 is some characteristic value that generally depends on 
B z , the instability is relatively weak and grows on the Alfven 
timescale characterized by the axial field. In contrast, if a > a c , 
the instability is much more efficient, and the growth time is of 
the order of the Alfven timescale for the toroidal field. 

In Fig. 8, the critical value of a c above which the system 
is unstable is plotted as a function of s for a range of verti- 
cal wavevectors relevant to stellar conditions. In fact, the local 
approximation in the axial direction applies if \k,s\ > 1 that is 
equivalent \q\ > 1 . On the other hand, to neglect gravity can be 




Fig. 8. The critical value of a that determines the onset of insta- 
bility as a function of e for \q\ = 10. Solid and dashed curves 
correspond to q — 10 and q = -10, respectively. The dash-and- 
dotted curve show critical a for m = 0, which does not depend 
on the sign of q. 



justified if \q\ > N/u>bo (see Eq. (11)), which generally imposes 
a stronger restriction on q since the ratio N/ojbq is typically > 1 
in stars. Therefore, in determining a c , we assume that the al- 
lowed \q\ should be large enough and we choose \q\ > 10, which 
corresponds to axial wavelengths shorter than A C /\Q. Since the 
presence of an axial field breaks symmetry between positive and 
negative q, critical curves are different for positive and nega- 
tive wavevectors. The region of a above the lines corresponds to 
configurations that are unstable for a given m. It turns out that 
true a c , which determines instability, corresponds to perturba- 
tions with positive m and negative q as shown in Fig. 8 by solid 
lines. Critical a decreases with decreasing s everywhere except 
a region of small s where the dependence a c (s) is very sharp. 
This particularly concerns the curve m = 0, which goes up very 
sharply at s < 0.01 and reaches the value 1 at s = in agreement 
with the result by Tayler (1973). However, the scale of Fig. 8 
does not let us see this sharp behaviour. Unfortunately, our code 
does not allow to follow the behaviour of critical a when it ap- 
proaches the value - 1 because of the singular character of the 
last term on the left hand side of Eq. (15). The corresponding re- 
gion is marked by crosses in Fig. 8. A more refined consideration 
is needed for this case which should perhaps include dissipation 



4. Discussion 

We have considered the linear stability of magnetic configu- 
rations containing the toroidal and axial fields, assuming that 
the behavior of small perturbations is governed by equations 
of the non-dissipative incompressible magnetohydrodynamics. 
This approximation is justified if the magnetic field is subther- 
mal and the Alfven velocity is low compared to the sound speed. 
The stability of magnetic configurations is a key issue for under- 
standing the properties of various astrophysical bodies such as 



peculiar A and B stars, magnetic white dwarfs and neutron stars. 
The magnetic instability can alter qualitatively the properties of 
configurations generated, for example, by dynamo in stars. Many 
dynamo models predict that the toroidal field should typically be 
stronger than the poloidal one, but such configurations can be 
unstable if the generated toroidal field does not decrease enough 
rapidly with s. The instability generates large- and small-scale 
motions that should alter the geometry of a generated magnetic 
field. 

Even though the poloidal field is weaker than the toroidal one 
in a number of dynamo models, its effect on the stability proper- 
ties cannot be neglected. This particularly concerns the behavior 
of the nonaxysimmetric perturbations considered in the present 
paper. If B z is relatively weak (B z < B^) then, typically, there 
exists a wide range of the azimuthal wave-numbers m for which 
the instability may occur. For any given m, only perturbations 
within some particular range of the vertical wave-vectors k z can 
be unstable. The growth rate is maximal for perturbations with 
k z of about 

m Bu, 



s B 7 



(19) 



Equation (19) is approximately equivalent to the condition that 
the Alfven frequency is vanishing somewhere within the cylin- 
drical layer (see Eq. (17)). It should be noted that the ratio B z /B^ 
can be rather low in stars, and the maximum growth rate at given 
m corresponds to very short axial wavelengths, 



2tt 



m B w 



(20) 



Taking, for example, B z /B v> ~ 10 and assuming that s is com- 
parable to the stellar radius, s ~ 10 11 cm, we find that the 
most rapidly growing modes should have the axial wavelength 
~ 6 x 10 9 /m cm. From our results it follows that the maximum 
growth rate increases slowly with increasing m and, therefore, 
perturbations with a very short azimuthal wavelength (very large 
m) should dominate the development of instability. For instance, 
the most rapidly growing mode with m — 100 corresponds to the 
axial wavelength ~ 6 x 10 7 cm, which is very short compared 
to the radius. Therefore, the instability of magnetic configura- 
tions can often be determined by the modes with very large m 
and an extremely short wavelength in the z-direction. This fact 
can cause problems in numerical modeling of the instability be- 
cause simulations will require a very high resolution in the <p- 
and z-directions. 

Depending on the profile of the toroidal field and the strength 
of the axial field, the instability can arise in two essentially dif- 
ferent regimes. In the case of a weak axial field, B^ » B z , the 
value of a that distinguishes between the regimes is « -1/2. If 
a > -1/2, then the instability grows on the Alfven timescale de- 
termined by the toroidal field and is rather fast. In this case, the 
growth time is 



r ~ OAp^suB, 



<p3 



yrs, 



(21) 



where su = s/10 n cm. If a < -1/2, then the growth time is 
given by same expression (18) but where B^ should be replaced 
by B z j = B z /10 3 G. Since B v » B z , the instability is slower and 
grows on the timescale determined by the axial field in this case. 
The transition between two regimes occurs at larger a if the axial 
field increases. 

It is fairly difficult to compare our results obtained for a sim- 
ple model with the available numerical simulations, which usu- 
ally use completely different basic magnetic configurations. For 



example, in calculations by Braithwaite (2006, 2007), the basic 
configuration was assumed to be either purely toroidal or purely 
poloidal, and the stability properties of such configurations dif- 
fer qualitatively from those considered in this paper. Recently, 
Braithwaite & Nordlund (2006) and Braithwaite (2008) have 
considered stability of the magnetic configuration with random 
initial fields. A vector potential was set up as a random field 
containing spatial scales up to a certain value. This random field 
was then multiplied by some screening function, so that the field 
strength in the atmosphere was negligible. 

This initial configuration contains both the toroidal and 
poloidal fields but is very different from our simple model. 
Nevertheless, some features seem to be in common even for 
such different models. Braithwaite & Nordlund (2006) and 
Braithwaite (2008) find that the instability can lead to different 
equilibrium configurations depending on the screening function. 
If the screening function for random fields decreases slowly or 
does not decrease at all, then the final equilibrium magnetic con- 
figuration is essentially non-axisymmetric. In contrast, an equi- 
librium configuration is closer to axisymmetry (but not axisym- 
metric) if the screening function decreases rapidly. This depen- 
dence on equilibrium configurations obtained in numerical cal- 
culations can reflect the regimes of "strong" and "weak" instabil- 
ities that correspond to different growth rates depending on the 
value of a. In accordance with our analysis, the growth rate at 
given e is higher for higher a (compare, e.g., Figs. 3 and 7). Our 
parameter a mimics to some extent the screening parameter p 
introduced by Braithwaite (2008) with decreasing a correspond- 
ing to increasing p. Therefore, we can expect from our analysis 
that non-aximetric instabilities should be more efficient for the 
screening function with p — than with p — 1 and that the final 
configuration exhibits stronger departures from axisymmetry for 
smaller p. This conclusion seems to be in qualitative agreement 
with the results of Braithwaite (2008). 

Our simple model does not take into account the stratifica- 
tion that can be important in many astrophysical applications. 
Basically, stratification provides a stabilizing influence if the 
temperature gradient is sub-adiabatic. However, this influence 
is small if perturbations have a relatively short wavelength in 
the axial direction, A < A c , such that inequality (11) is satisfied. 
Our results are related to this case. The case when A > A c and 
stratification is important will be considered elsewhere. 
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